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We analyze the behavior of a suspension of rigid rod-like particles in shear flow using a mod- 
ified version of the Doi model, and construct diagrams for phase coexistence under conditions of 
constant imposed stress and constant imposed strain rate, among paranematic, flow-aligning ne- 
matic, and log-rolling nematic states. We calculate the effective constitutive relations that would be 
measured through the regime of phase separation into shear bands. We calculate phase coexistence 
by examining the stability of interfacial steady states and find a wide range of possible "phase" 
behaviors. 
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I. INTRODUCTION 

Shear flow has profound effects on complex fluids. 
It can perturb equilibrium phase transitions, such as 
the isotropic-to-nematic (I-N) liquid crystalline tran- 



thcrmotropic melts 
the ncmatic-smectic 



sition in wormlike micelles 
[Q-Q, or rigid-rod suspensions"^ 
transition in thermotropic liquid crystals ]10fl; and the 
isotropic-to-lamellar transition |p"l|| in surfactant sys- 
tems. Shear can also induce structures, such as the well- 
known multi-lamellar vesicles (onions) in surfactant sys- 
tems p^-[j^], that exist only as metastable equilibrium 
phases. Another well-known effect is the transition be- 
tween orientations of diblock copolymer lamellae in either 
the steady shear flow [jl5 16 1, or the oscillatory shear flow 
0-|l9| , as a function of shear rate or frequency, and tem- 
perature. 
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FIG. 1. Stress-strain-rate curves for the Doi model with 
different excluded volume parameters u (taken from Fig. |^ 
below). The dashed line segments are unstable (unphysical) 
steady states. The straight lines indicate possible coexistence 
between states I and II under conditions of common stress 
(horizontal lines) or strain rate (vertical line). 

A related phenomenon is dynamic instability in non- 
Newtonian fluids whose theoretical homogeneous stress- 



strain-rate constitutive relations exhibit multi- valued be- 
havior, as in theories of polymer melts 1 20 2l[] and worm- 
like micelles [^2|-f24f . Such models may describe, for ex- 
ample, the spurt effect, whereby the flow rate of a fluid 
in a pipe changes discontinuously as a function of applied 
pressure drop p5| . A non-monotonic constitutive curve 
as in Fig. |l| typically has a segment (shown as a broken 
line) where bulk flow is unstable. If a mean strain rate 
is imposed which forces the system to lie on an unstable 
part of the constitutive relation, a natural resolution for 
this instability is to break the system into two regions, 
often called bands, one on the high strain rate branch and 
one on the low strain rate branch, to maintain the overall 
applied strain rate. The most important unresolved ques- 
tion about these banded flows is, what determines the 
stress at which the system phase separates into bands? 
Experiments on many systems (reviewed in Sec. HA), 
particularly the wormlike micelle surfactant systems, re- 
veal that there is a well-defined and reproducible selected 
stress in a wide class of systems 

There have been many suggestions for determining the 
selected stress. Some workers have assumed the existence 
of a non-equilibrium potential and a variational principle 
p6| , p7| . This possibility is intriguing, although it remains 
unproven. Early studies postulated a jump at the top of 
the stable viscous branch ("top jumping") [2l]j2^ ,|2^| , but 
experiments have shown that this is not the case p9]. Re- 
cent studies have solved the homogeneous flow equations 
in various geometries using sophisticated hydrodynamic 
flow-solvers and found a selected stress |^,^l[. How- 
ever, evidence is growing |3^] that these calculations have 
history-dependent stress selection (which is in fact no se- 
lection) or introduce gradient terms due to the discretiza- 
tion of the system. A final method, which we follow here, 
has been to incorporate (physically present) non-local 
contributions to the stress 32 3^ [57| , and ex- 



amine the equations of motion under steady banded flow 
conditions. 

Here we extend previous work ||] and calculate phase 
diagrams for rigid-rod suspensions in shear flow, solving 
for the interfacial profile between phases and using its 
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properties to determine the coexistence stress. As Fig. [T] 
indicates, phase separation is possible at either a speci- 
fied stress (horizontal tie lines) or a specified strain rate 
(vertical tie lines). Only recently has the latter possibil- 
ity been speculated upon 1 28 |,|27j , and found experimen- 
tally [p8[ . We explore this possibility explicitly for our 
model system, which possesses, in addition to the high 
and low strain rate (paranematic and nematic, respec- 
tively) branches shown in Fig. [I], a second high strain rate 
branch in which the rods stand up in the flow, parallel to 
the vorticity direction, instead of lying in the shear plane 
p9| . We study coexistence with this so-called 'log-rolling' 
phase and find a rich non-equilibrium phase diagram. 

The summary of this paper is as follows. In Section II 
we discuss the general issues of shear banding and phase 
separation in flow, and summarize the primary exper- 
imental evidence for this behavior. In Section III wc 
present the modified Doi model ^(],|d|] and in Section IV 
we briefly discuss our algorithm for calculating the phase 
diagram. The general aspects of the interface construc- 
tion will be discussed elsewhere 35 . We present the 



results for common stress and strain rate phase separa- 
tion in Section V and VI, respectively, and discuss some 
of the implications for metastability and experiments un- 
der controlled stress or controlled strain rate conditions. 
We finish in Section VII with a discussion and summary. 
While some of these results have been briefly summa- 
rized elsewhere [Q , the current paper is a complete and 
self-contained discussion of the problem. 

The reader interested in the phenomenology of phase 
diagrams for sheared complex fluids rather than liquid 
crystals may safely skip Section III; the rest of the pa- 
per is general, and much of the discussion applies to any 
system undergoing phase separation in shear flow. There 
are, essentially, two steps to calculating phase behavior 
in flow. One must derive the dynamical equations of mo- 
tion for fluid flow, composition, and the relevant struc- 
tural order parameter(s), which is quite difficult. Then, 
one must understand how to solve them and interpret 
the results. While the modified Doi model does not ex- 
haust all possible phase diagrams (in particular, a shear- 
thickening model would be a nice complement), it has 
many universal features. One extremely important con- 
cept is that density and field variables are ill-defined in 
non-equilibrium systems: either stress or strain rate may 
act as a control parameter analogous to an equilibrium 
field variable (e.g. pressure, chemical potential), corre- 
sponding to the different orientations of the interface be- 
tween coexisting phases. Also, one can gain much intu- 
ition from the underlying stress-strain-rate-composition 
surface, a fact which we feel has been underappreciated 
until now. 



II. SHEAR BANDING 



A. Experimental Evidence 



Shear banding has been confirmed in many systems 
through direct optical and NMR visualization, and de- 
duced from rheological measurements. The best-studied 
systems are surfactant solutions of various kinds, includ- 
ing wormlike micelles and onion-lamellar phases. Rehage 
and Hoffmann pi| measured a plateau in the stress- 
strain-rate relation for wormlike micelles in shear flow. 
This behavior has since been seen in a number of worm- 
like micellar systems in various flow geometries, by the 
Montpellier (ML Strasbourg g|, Edinburgh ||, and 
Massey groups |44|-|47[| . Berret et al. || visualized shear 
bands in the plateau region of the stress-strain-rate 
curves using optical techniques, providing proof of band- 
ing; and Callaghan et al. ]4^-|47|] used NMR to measure 
the velocity profile in various geometries (including Cou- 
ette, cone-and-plate, and pipe geometries). 

The transition in these cases is to a strongly-aligned, 
possibly nematic, phase of wormlike micelles which has a 
lower viscosity than the quiescent phase. It is not known 
how the length distribution changes in flow, although 
this is certainly an important aspect of these 'living' 
Wormlike micellar system can possess an 



systems 

equilibrium nematic phase, and in some cases the shear- 
induced phase is obviously influenced by the proxi mity 
of an underlying nematic phase transition [[ip|4^,|3]j30f . 
However, many wormlike micellar systems undergo band- 
ing at compositions much more dilute than that for I-N 
coexistence, and it is probable that in these cases flow 
instability is due to the non-linear rheology of these sys- 
tems, which is in many respects similar to that of the 
Doi-Edwards model of polymer melts p^] . Since there 
are at lease two possible effects (a nematic phase transi- 
tion and flow-instability of the micellar constitutive re- 
lation) apparently leading to flow-instability, these sys- 
tems are quite rich. It is tempting to analyze the extent 
to which these systems display behavior analogous to the 
kinetics of equilibrium phase separation, and groups have 
recently begun to study the kinetics of non-equilibrium 
phase separation P, p9|]5l] ]. 

Pine and co-workers have recently studied a wormlike 
surfactant system at extreme dilutions and found, sur- 
prisingly, that for low enough concentrations (but still 
above the overlap concentration) shear induces a vis- 
coelastic phase that they interpret as a gel |52]-Q. The 
origins and structure of this gel are currently unknown. 
In controlled stress experiments they observe shear band- 
ing and a 'plateau' for stresses higher than a certain 
stress, in which the strain rate decreases as shear induces 
the gel. Above the stress at which the gel fills the sample 
cell the strain rate increases again to complete a dramatic 
S curve. For controlled strain rate experiments the sys- 
tem jumps, at a well defined strain rate, between the gel 
and solution phases. 

Another well-studied system is the onion lamellar sur- 
factant phase, originally studied by Roux, Diat, and Nal- 
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let |lJ-|TJ]. These systems display a bewildering variety 
of transitions between lamellar, aligned-lamellar, onion, 
and onion crystal phases of various symmetries, as func- 
tions of applied shear flow, temperature, and composi- 
tion. As an example, one particular system undergoes 
transitions, with increasing strain rate, from disordered 
lamellae to onion, to onion-lamellae coexistence (in which 
coexistence is inferred from a plateau in the stress-strain- 
rate curves), to well-ordered lamellae {f|. Recently Bonn 
and co-workers |38| found shear-induced transitions be- 
tween different gel states of lamellar onion solutions with 
shear bands (visualized by inserting tracer particles) ori- 
ented with interface normals in the vorticity direction, 
indicating phase separation at common strain rate in- 
stead of common stress, as we clarify below. In this case 
the averaged stress-strain-rate constitutive relation fol- 
lowed a sideways S curve under controlled strain rate 
conditions. 

Mather et al. j?]] have recently studied a thermotropic 
polymer liquid crystal using visual and rheological mea- 
surements, and inferred a shear-induced nematic phase 
transition and phase separation, the latter which they 
attribute to polydispersity. 

In summary, shear-banding has been seen in several 
systems, and in all cases is associated with some flow- 
induced change in the fluid microstructure. Most sys- 
tems are still poorly-understood (5^JTJ,Q and, given the 
range of complexity, it is certain that many qualitatively 
new phenomena remain to be discovered. 



B. Theoretical Issues 

The crux of the problem from a theoretical point of 
view may be appreciated from Fig. [I]. These stress- 
strain-rate curves are somewhat reminiscent of pressure- 
density (p — p) isotherms for a liquid-gas system. Curve 
segments with negative slope, da xy /dj < 0, are unstable 
and cannot describe a physical state of a bulk homo- 
geneous system. Analogously, isotherms with negative 
slopes dp/ dp < have negative bulk moduli and are 
unstable. The liquid-gas system resolves this instabil- 
ity by phase-separating into regions of different densi- 
ties (according to the lever rule to maintain the average 
density). Similarly, the banded flows seen in the exper- 
iments described above appear to be a non-equilibrium 
phase separation into regions of high and low strain-rate, 
maintaining the applied mean strain rate. 

In previous work J^,^,§| we constructed a 'phase dia- 
gram' by pursuing an analogy between homogeneous sta- 
ble steady states and equilibrium phases. As in equilib- 
rium, non-equilibrium 'phases' may be separated, in field 
variable space, by hyper-surfaces representing continu- 
ous (e.g. critical points/lines) or discontinuous ('first- 
order') transitions. Coexistence implies an inhomoge- 
neous state spanning separate branches of the homo- 
geneous flow curves. Note, however, that there is an 



ambiguity in connecting separate branches of the homo- 
geneous flow curves in Fig. [j]. The top curve permits 
coexistence of states with the same stress and different 
strain rates, while the lower curve also allows coexistence 
of states with the same strain rate and different stresses 



Fig. U shows that phase separation at a common stress 
occurs such that the interface between bands is parallel 
to the vorticity- velocity plane (annular bands, in Couette 
flow), while phase separation at a common strain rate 
occurs with the interface between bands parallel to the 
velocity-velocity-gradient plane (stacked discs, in Cou- 
ette flow). 





FIG. 2. Geometries for phase separation at common stress 
(left) or strain-rate (right) in a Couette rheometer. For phase 
separation at a common stress (left) phases / and II have 
different strain rates, while at a common strain rate (right) 
they have different stresses, z is the vorticity axis, x is the 
flow direction, and y is the flow gradient axis. 

This highlights a striking contrast between equilibrium 
and non-equilibrium systems. In equilibrium the field 
variables (pressure, temperature, chemical potential) are 
uniquely defined and determine phase coexistence. In 
sheared fluids, one needs an extra field variable to deter- 
mine the extended phase diagram. However, for a system 
with more than one choice of coexisting geometry, the ap- 
propriate held variable may not necessarily be identihed 
a priori. The complete answer of how to determine (the- 
oretically) the dynamic field variable is not known. Of 
course the nature of the constitutive relation may help, 
for example the top curve of the Fig. [jjdoes not allow the 
strain rate as the field variable. We will come back to dis- 
cuss some possib le answers to this interesting problem in 
subsection VII B (see also |^8| for other suggestions). 

Another important difference from equilibrium sys- 
tems is evident when, say assuming the system choose 
to form shear bands at common stress, we try to deter- 
mine at which stress a system shear bands. The con- 
stitutive relations shown in Fig. [I] are calculated for ho- 
mogeneous states, and there is no apparent prescription 
for determining the selected banding stress, despite the 
experimental evidence for a selected stress. A similar ap- 
parent degeneracy occurs in first order phase transitions 
in equilibrium statistical mechanics, but is easily resolved 
by demanding that the system minimize its total free en- 
ergy, or, equivalently, by appealing to the convexity of 
the free energy of the equilibrium thermodynamic sys- 
tems |55j . This leads to equality of field variables between 
two phases and the common tangent condition (e.g. the 
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Maxwell equal areas construction for liquid-gas coexis- 
tence |5(| , or the equal osmotic pressure condition, aided 
by equal chemical potential, in rod suspensions J57| ). 

In the shear band problem, an unambiguous resolution 
of this degeneracy is to consider the full inhomogeneous 
(i.e. non-local in space) equations of motion, and deter- 
mine phase coexistence by that choice of field variables 
(appropriately chosen by hand) for which there exists a 
stationary interfacial solution to the steady-state differ- 
ential equations of motion For zero stress this 
technique reduces, as it should, to minimization of the 
free energy. The importance of inhomogeneous terms 
in fluid equations of motion has been noted by several 
groups, who pointed out that the standard fluid equa- 
tions can have ill-defined mathematical solutions 58 if 
such terms are not included. Of course, if the phase di- 
agram depends sensitively on the form or magnitude of 
the inhomogeneous terms one need a detail understand- 
ing of the underlying physics. The use of a stable inter- 
face to select among possible coexisting states was first 
postulated for non-linear dynamical systems, as far as 
we know, by Kramer fl59[ , and later by Pomeau pCfl ; and 
first applied (independently) to complex fluids in Ref. Q . 
The inclusion of gradient terms in constitutive relations 
is rapidly gaining acceptance, as recent preprints by Gov- 
eas (phase separation of model blends of long and short 
polymers) |p6[ and Dhont (introduction of model gradi- 
ent terms to resolve stress selection) Jl?} indicate. 

In this work we study a model for rigid-rod suspensions 
in shear flow. While there are certainly ongoing exper- 
iments on these systems Q|, the primary motivation for 
this extended work is to explore the manner in which 
phase separation and coexistence occurs in complex flu- 
ids in flow. The approximations used in obtaining our 
equations are severe (including a decoupling approxima- 
tion whose defects are well-known J6l| ), and we expect 
qualitative agreement at best. However, this is the first 
complete study of which we are aware of non-equilibrium 
phase separation of a complex fluid in flow for a concrete 
model, and we hope it illuminates the phenomenology of 
flow-induced phase transitions. 



III. METHODOLOGY 

We seek the equations of motion for a solution of rod- 
like particles. The most useful dynamic variables describ- 
ing the long-wavelength hydrodynamic degrees of free- 
dom are the volume fraction 0(r), the fluid velocity v(r), 
and the nematic order parameter tensor 

Qa/3(r) = {VaVp - \5 a p), (3.1) 

where u is the rod orientations, (•) denotes an average 
around the point r. Previous studies of liquid crystals 
under shear flow have been either for thermotropics , 
where the issues we present below associated with compo- 
sition coupling are not present; or homogeneous suspen- 
sions f§j39], where phase coexisting was not considered. 



Our work below is based on the model extending that 
of Doi (HQ. See, et al. g studied the Doi model in 
shear flow, but did not attempt to consider phase coex- 
istence. Bhave, et al |]39[] analyzed this model in more 
detail, but did not consider realistic phase separation 
behavior. We augment this model with reasonable es- 
timates for translational entropy loss upon phase sepa- 
ration and for the free energy cost due to spatial inho- 
mogeneities. Zubarev studied shear-induced phase sepa- 
ration in a variation of the Doi model in flow based on 
the equality of non-equilibrium free energies, calculated 
from the flow-perturbed orientational distribution func- 
tion pq |. Zubarev only considered phase separation at 
a common strain rate, and did not treat the rheological 
response (stress) of the system or log-rolling states. 



A. Equations of Motion 

The free energy (e.g. as in Ref. |5^]) is given by 

(1- 



F(4>, Q)= k B T / cPr 



log< 



1ok(1- 



\ (1 - \u) Tr Q 2 - ±uTr Q 3 + \u (Tr Q 2 )' 



If <v*) a 



(3.2) 



Here, "Tr" denotes the trace, v r and v s are rod and sol- 
vent monomer volumes and 



u — u 2 cdL 2 1 



(3.3) 



is Doi's excluded volume parameter |lO],[n| , where c is the 
concentration (number/ volume) of rods of length Lq and 
diameter d, and v 2 is a geometrical prefactor (Ref. ^] 
estimated v 2 — 57r/16 ~ 0.98). The volume fraction (f> is 



in terms of which 



4> = cv r , 



u = <pL — , 
a 



(3.4) 



(3.5) 



where L = Lo/d is the rod aspect ratio and a is an 0(1) 
prefactor defined by 



ad Lq. 



(3.6) 



For spherocylinders, a = tt[1 — l/(3L)}/8 which reduces, 
in the limit L — > oo, to a = ir/8 ~ 0.39. We use u and 4> 
interchangeably below as a composition variable. 

In much of what follows we make two further assump- 
tions to reduce the number of parameters in our model. 
We fix v s by assuming 



v r = Lv* 



(3.7) 
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which corresponds to a particular volume of the solvent 
molecules relative to that of the rod-like molecules. Fur- 
ther, we assume that the geometric factor v-i /a has the 
value unity, so that 



u = <pL, 



(3.8) 



which corresponds to a particular shape of the rigid-rod 
molecules. These two assumptions specify the detailed 
shape and volume ratio of the system we study below. 
For slightly different systems with v r Lv s 1, 
our work should still provide an accurate qualitative pic- 
ture. 



The first two terms of Eq. ( ^2 ) comprise the entropy 
of mixing, and the first three terms in square brackets 
are from Doi's expansion of the free energy (derived per 
solute molecule) in powers of the nematic order parame- 
ter Q. These terms were derived from the Smoluchowski 
equation for the distribution function of rod orientations 
We keep the expansion to fourth order to de- 
scribe a first order transition and give the correct quali- 
tative trends. 

Assuming Eqs. (3.7-3.8), we calculate the following 



biphasic coexistence regions, 

{m = 2.6925, mat = 2.7080} 
{m = 2.6930, mat = 2.7074} 



(L = 5.0) (3.9) 
(L = 4.7), (3.10) 



where uj and un are the excluded volume parameters 
(compositions) for the coexisting isotropic and nematic 
phases, respectively. Note the very weak dependence of 
the biphasic regime (in the sc aled variable u = L<f>) on L. 

The last two terms in Eq. (3^) penalize spatial inho- 
mogeneities. By adding the single term proportional to 
K we have assumed a particular relation for the Frank 
constants, {K x =K 2 = K, K 3 = 0) . Although Odijk 

has calculated these constants for model liquid crystals 
(in the nematic regime) |34j, we will see below that this 
choice is probably unimportant for this model. More gen- 
erally, we expect the Frank constants to vary as functions 
of Q(r) in physical systems, a situation which we have 
not addressed here. The final term penalizes composition 
gradients [[35| . We are not aware of any calculations of g 
for solutions of rod- like particles. In Eq. (3.2), we assume 
an athermal solution with no explicit interaction energy. 

The nematic order parameter obeys the following equa- 
tion of motion ■ 



(S t +v-V)Q = F( K ,Q) + G(0,Q) 



(3.11) 



where K a p = Vpv a . In Eq. ( |3.1l| ) the (reactive) ordering 
term F is given by 

F{h,Q)^Ik s + k i -Q + Q-h t -2(Q + \I)Tt{Q-k) 1 

(3.12) 

where k s is the symmetric part of k and J is the identity 
tensor. For simplicity, we have chosen the form appro- 
priate for an infinite aspect ratio (the prefactors differ 



by 0(1) constants for finite aspect ratios jyj). The cou- 
pling F to the flow both induces order, and dictates a 
preferred orientation. The dissipative portion G is 



where 



D r = 



(l-|TrQ 2 )2( c L3)2- 



(3.13) 



(3.14) 



is the collective rotational diffusion coefficient and 



H 



5T_ 
SQ 



8T 



- ITr 

3 SQ 



(3.15) 



is the molecular field. D ro is the single-rod rotational 
diffusion coefficient and 1^1 is an 0(1) geometrical pref- 
actor, which will be fixed below Eq. 3.30. The rotational 
diffusion coefficient is 



D r — 



k B T\nL 



(3.16) 



where r\ is the solvent viscos ity. The Q-dependence in 
the denominator of E q. (3.14 ) enhances reorientation for 
well-ordered systems |4C |. Our choice for D r is crude, 
since it applies to rods in concentrated solution and 
we use it in the concentrated and semi-dilute regimes. 
As with many of our approximations, this gives us a 
tractable model system with which to study the phe- 
nomenology of phase separation. 



Doi and co-workers derived Eq. (3.11) for homogeneous 
systems. We extend this to inhomogeneous systems by 
including the gradient terms implicit in the functional 
derivative which defines H. Our choice of F is the 
so-called quadratic closure approximation to the Smolu- 
chowski equation This approximation ensures that 
the magnitude of the order parameter remain in the phys- 
ical range in the limit of strong ordering, but is known 
to incorrectly predict phenomena such as director tum- 
bling and wagging. Many workers have investigated the 
subtleties of various closure approximations and the de- 
gree to which they reproduce realistic flow behavior |n) . 
Since our primary goal is to explore the method for cal- 
culating phase behavior and outline some of the possibil- 
ities for coexistence under flow, we confine ourselves to 
this well-studied model. 

The fluid velocity obeys E3j4l]J66| : 



p (d t + v • V) v = V ■ [2 V k s + tr{4>, k, Q)] 



5F, 

6(f) 



- Vp, 
(3.17) 



where rj is the solvent viscosity, p the fluid mass density, 
and the pressure p enforces incompressibility, V-v = 0. 
For the low Reynold's number situations considered here, 
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and for steady shear flow, we will equate the left-hand 
side of the equation above to zero. 

The constitutive relation for the stess tensor <r(</>, k, Q) 
was derived by Doi and co-workers, and includes dissipa- 
tive and elastic parts. Since the elastic stress dominates 
pOf, we keep only this part: 



& — & elastic 

= 3H + H Q Q H -VQ af3 



ST 



(3.18) 



The first term of Eq. ( J3 . 1 8| ) was given by Doi [|40| , while 
the last three terms were derived later Q and are equiv- 
alent to the elastic stress due to Frank elasticity J63[ , 
generalized to a description in terms of the nematic or- 
der parameter Q rather than the nematic director. Note 
that the last three terms vanish for a homogeneous sys- 
tem. 

Finally, the composition equation of motion is of the 
Cahn-Hilliard form §1; 



(d t + v • V) <j> = -V • J 

= V-M-V/i, 



(3.19) 



where M is the mobility tensor and the chemical poten- 
tial is given by 



ST 



V = 



(3.20) 



The diffusive current is J = — iW- V/i. The compl ete d y- 
namics is thus described by Eqs. (3.11,3.17) and (3.19). 

The dynamical equations of motion for other com- 
plex fluids have the same theoretical structure: equa- 
tions of motion for the conserved quantities and the 
broken-symmetry or flow-induced structural order pa- 
rameter (analogous to Q), and a constitutive relation for 
the stress as a function of composition and order pa- 
rameter f36| . For a given system and set of equations 
of motion, the analysis below is gene ric. For some local 
models, internal dynamics (Eq. |3.1l| ) can be eliminated 
to give the stress as a history integral over the strain 
rate. In polymer melts |2(J , and in wormlike micelles j2j| 
far from a nematic regime, this leads to non-monotonic 
stress-strain-rate curves. However, augmenting these in- 
tegral theories with non-local terms to calculate interface 
profiles is non-trivial. 



B. Steady-state conditions 

In this work we study planar shear flow, specified by 
dv x (r) 



dy 



= 7(r). 



(3.21) 



For homogeneous flows v(r) = 7yx. The phase diagram 
is given by the domains of stable steady-state solutions 



to the equations of motion for applied shear stress or 
strain-rate, in the phase space spanned by 

(4>,<T xy ), (common stress) (3.22) 

(0,7), (common strain rate). (3.23) 

For phase separation at common stress the stress is uni- 
form and the strain rate partitions between the two 
phases; while for phase separation at common strain rate 
the strain rate is uniform and the shear stress partitions 
between the two phases. 

The strain rate tensor is given by 




Upon rescaling, 



7 L 2 



&D ro Vil>2 



3k B T 



(3.24) 

(3.25) 
(3.26) 



the stead y-state condition for the order parameter 
(Eq. |TL|) is 



= 



u 2 L 2 (l - §TrQ 2 ) 2 



(3.27) 



where F = jF and 



H 



1--J Q-u \Q --TtQ 



uQTtQ 2 - KV 2 Q. 

(3.28) 



In steady state planar shear flow the velocity gradi- 
ents are normal to the flow direction, so the convective 
derivative vanishes and Eq. (3.27) specifies the order pa- 
rameter in a homogeneous flow. Under the se co nditions 
integration of the momentum equation Eq. ( (3 .17 ) gives a 
constant stress, 



<t = er - pi + 2l]K S , 



(3.29) 



where ctq is the boundary stress. The rescaled shear 
stress is 



Aj — uL 



H + K(V 2 QQ-QV 2 Q) 



xy 

(3.30) 



where A = 2zyi^|(lnL)/(37r) is a constant of order unity: 
we take A = 1 for the remainder of this work, which 
corresponds to a particular choice for v\. As with the 
assumpt ions of molecular geometry embodied in v 2 and 
a (Eqs. 3.7-3.8), different values for A should not quali- 



tatively change the nature of our results. 

Integrating the steady-state composition equation 
fl3~19l ) and using the boundary condition that material 
cannot enter or leave the system, we find 



G 



Mo = MO) 



(3.31) 
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L(l 



In (1-0)] 



+(j)L—F Doi + \K (VQ) 2 - gL\7 2 cj) (3.32) 



where Eqs. ( j3~7l - |3.8| ) have been used to specify the molec- 
ular geometry, /io is a constant of integration, and 



F D oi = \ (1 - \v) Tr Q 2 - \uTi Q 3 + \u (TtQ 



(3.33) 



Note that the mobility tensor M plays no role in the 
steady-state conditions, or in the resulting phase dia- 
gram. 

Eqs. ( 3.27 3.30| ) and ( 3.31 ) completely specify the sys- 
tem in planar shear. Solving these equations will occupy 
the remainder of this work. Note that variables 7, <r, 
fi/k B T are all dimensionless quantities. 



IV. CALCULATION OF PHASE DIAGRAMS 



A. Interface Calculation 



The phase diagram is specified by solving 
Eqs. ( |3~2^ , |3~30| , p3ll ) for given fiQ and boundary stress 



Non-equilibrium 'phases' are defined as the stable 
steady-state space-uniform solutions to these equations. 
These inhomogeneous equations comprise a set of ordi- 
nary differential equations, through the gradients that 
appear in the stress and in the functional derivatives that 
define fj, and H. The only parameters of the theory are 
the rod aspect ratio L and the ratio of elastic constants, 
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FIG. 3. Homogeneous stress a xy vs. strain rate 7 behavior 
for various excluded volumes, L — 5.0 and A = 1.0. Dotted 
lines mark unstable branches. Similarly, curves along which 
d/j,/d(j> < are linearly unstable. 

We first fix cf> ( u) an d solv e the homogeneous algebraic 
versions of Eqs. ( 3.27 ) and ( 3.30 ) for Q and 7 as a func- 
tion of cr^y* ■ This is done for all <p. Because !F((j>, Q) de- 
scribes an I-N transition, at a given stress, multiple roots 
exist, with distinct strain rates and Q. Fig. |^ shows the 
stress stra in-rate rela tions for homogeneous solutions to 
Eqs. ( |3.27|) and (|3.30| ) for L = 5.0 and A = 1.0. 




FIG 4. Stress-strain-composition surface for the curves in 
Fig |. The plane is at d° xy = 0.05 

The isotropic branch has a larger viscosity than the 
nematic branch, and has an increasing effective viscosity 
for increasing concentration, reflecting the contribution 



(K may be absorbed into the length scale of the system). 



uH in Eq. (3.30). Conversely, the nematic branch has a 
lower stress at higher concentrations due to the increased 
nematic order which permits less-hindered motion. 

For a dilute isotropic system (curve a), shear flow con- 
tinuosly induces nematic order. A more aligned system 
has a lower effective viscosity, so the stress a(j) increases 



*In the few cases where the phase diagram in the a xy ~fi 
plane has a transition line parallel to the fi axis, one must 
first fix /io, and then determine ctq. 



7 



slower than linearly (shear-thins) as the magnitude of the 
order parameter Q increases. Eventually the system at- 
tains, smoothly, a high strain rate state with a much 
lower viscosity than in the limit of zero stress. For more 
concentrated systems (curves b, c and d) shear flow in- 
duces a transition to a nematic phase with lower viscosity, 
and (7(7) is non-monotonic. There is a region of stresses 
for which two stable strain rates exist, on either the ne- 
matic or isotropic branches of the constitutive curve. For 
compositions inside the biphasic regime (curve e) both 
nematic and isotropic branches exist in the limit of zero 
stress, with the isotropic branch losing stability at high 
enough stress. Finally (not shown) for highly concen- 
trated systems only nematic branches exist. 

As mentioned in Sec. IIB, we calculate the phase dia- 
gram by explicitly constructing the coexisting interfacial 
solution ||. In common stress coexistence, for example, 
the coexisting states have different strain rates and, gen- 
erally, different compositions. Hence, they connect the 
high and low strain rate branches of two different curves 
in Fig. |^. It is easiest to visualize this by considering 
the intersection of a plane at a given stress <j° with the 
surface a xy (j,u), as in Fig. 0. 




2.59 
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FIG. 5. (a) Reduced strain rate j(u) and (b) chemical po- 
tential aj(m) for the stress contour in Fig. ^ {ffxy = 0.05). The 
tie line is calculated using the interface construction. 

At a given stress, the strain rate varies with composi- 
tion as shown in Fig. ||a. At coexistence, the chemical 
potential n(r) must be constant through the interface, as 



dictated by Eq. (3.31). The functional form of the non- 
equilibrium chemical potential is known from Eq. ( 3.32 ), 
and depends on the strain rate through the dependence of 
the nematic order parameter on the strain rate in steady 
state. We plot /j,(u) in Fig. ||b. There is a continuum 
range of fi, which allow possible coexisting pairs of states. 
(Recall that u is proportional to the rod volume fraction 
0). 

We now impose the interface solvability condition as 
follows. For a given stress a a xyl we determine a specific 
coexistence chemical potent ial up. which a llows a stable 
interfacial solution to Eqs. ( |3.27 . 3.30 13. 31 ). I n pra ctice, 
we eliminat e j(r ) from Eq. (3.27) using Eq. (3.30), and 
solve Eqs. ( |3.27| ) and ( 3.31 ) for the interfacial profile, 
with boundary conditions (fixed Q and u) chosen by two 
points on the low and high strain rate branches of Fig. ^Jd 
with the same /io- We adjust /j,q until a stationary interfa- 
cial profile is found. This solvability criterion give sharp 
selection on /z, and in this way determine a tie line on the 
7 — u plane, Fig. ||a. By varying the stress we compute 
the entire phase diagram in the a xy — u and 7 — u planes. 

For phase separation at a common strain rate the 
construction is analogous. One slices a vertical plane 
through Fig. |^ at a given strain rate, constructs the curve 
n(u) along the intersection with the surface, and searches 
for a stationary interfacial solution. 

The interface calculations are carried out by discretiz- 
ing the system on a one-dimension m esh and, from 
smooth initial conditions, evolving Eqs. ( 3.27 , 3.30 , 3.31 ) 
forward using fictitious dynamics calculated with an im- 
plicit Crank-Nicholson scheme. Spatial variations are 
only allowed in the direction in which phase separation 
occurs, so we replace 



d 

— — common stress 
oy 



— — common strain rate, 
. az 



(4.2) 



where z is in the vorticity direction. We fix the values 
of u and Q at either side of the interface to lie on the 
high and low strain rate branches of Fig. ||, begin with 
smooth initial conditions, and let the system "evolve" 
towards steady state. An interface develops between the 
two phases, and moves to one boundary or the other. For 
a given stress, coexistence is determined by that chemi- 
cal potential fi for which a stationary interface lies in the 
interior of the system (in the limit of large system size) 
B. An analogous construction may be made by main- 
taining a fixed mean strain rate on the unstable part of 
a homogeneous curve, and then starting up the system 
and allowing it to select a stress and chemical potential. 
In either case the selected stress is that stress for which a 
stationary interfacial solution between the high and low 
strain rate branches exists. Such an interfacial solution 
is known in dynamical systems theory as a heteroclinic 
orbit p9| , |60| js|] , and further work will investigate this in 
more detail for simpler model systems ]35||}g , g7] . 
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We restrict the nematic order parameter to 

/ 91 ?3 \ 
Q = <?3 <& , (4.3) 

\ -(qi+q 2 ) J 

since all steady state solutions with non-zero elements 
Qxz or Q yz are unstable due to the symmetry of shear 
flow (39). In a similar calculation for thermotropic ne- 
matics in shear flow we have found that this restriction 
on Q reproduces the same selected stress as that obtained 
when keeping the full tensor ||. 

For planar shear flow and a wide class of equations of 
motion we have shown that, if a coexisting solution ex- 
ists, it occurs at discrete points in the parameter set [ {35| . 
For example, for a given stress cr° y , coexistence can oc- 
cur only at discrete values for /no, that is, along lines in 
the field variable space spanned by a xy —fj,. This is anal- 
ogous to equilibrium systems where, for example, phase 
transitions in a simple fluid occur along lines, rather than 
within regions, in the pressure-temperature plane. 

Note that a one-dimensional calculation does not de- 
termine the stability of the interfacial solution with re- 
spect to transverse undulations (capillary waves), which 
could be important in, particularly, the common stress 
geometry Q. 




FIG. 6. Constitutive relations for I, N, and L states, for 
L = 5.0 and two values for the excluded volume parameter. 



The I phase is stable at lower volume fractions and 
merges with the N phase at high strain rates. The L 
phase is stable only at higher volume fractions, and is 
destabilized at high enough strain rates. For low strain 
rates the stress of the L state is lower than that of the N 
state, which is lower than that of the I state (see Fig. 



B. Homogeneous Solutions 



The modified Doi model in the quadratic closure ap- 
proximation has three stable solutions in homogeneous 
planar shear flow. We refer the reader to Bhave, et al. 
for further details fl39(l. 



I Paranematic: The paranematic state I induced 
from a disordered equilibrium phase. The order pa- 
rameter Q is small and fairly biaxial, with major 
axis lying in the shear plane at an angle of almost 
7r/4 relative to the flow direction. 



Fig. |^ shows the regions of stability of the various 
states. The loop in Fig. ^a occurs for compositions such 
that the constitutive curve cr( A f) has the shape of curve 
b in Fig. ||. Similar phase-plane plots were calculated by 
Bhave et al. |39| and See et al. H. They did not con- 
sider the mixing entropy needed to generate a realistic 
nematic transition, however, and always generated solu- 
tions for a given strain rate instead of a given stress, (this 
explains the absence of a loop in their phase-plane plot 
7 — </>)■ Their plots (compare Fig. 5 of Ref. (Ml) corre- 
spond to truncating the loop in Fig. 0a. FigT^b has a 
similar, barely discernable, loop near the critical point, 
within which there are no stable states. This instabil- 
ity is due to the instability of the composition equation, 
Eq. ( 3.1S ) . In this region 



Flow- Aligning Nematic: The flow-aligning nematic 
state is much more strongly-aligned, has slight bi- 
axiality induced by the flow, and has the major axis 
of alignment in the flow plane at an angle of a few 
degrees relative to the flow direction. The I and N 
states have the same symmetry. 



<0, 



(4.4) 



L Log-Rolling Nematic: The log-rolling phase is also 
a well-aligned and almost uniaxial phase, but with 
major axis of alignment in the vorticity (z) direc- 
tion, so the rods spin about their major axes. 



which is equivalent to a negative diffusion coefficient, and 
is analogous to the conventional definition of the spinodal 
line for ordinary equilibrium demixing. 
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FIG. 7. Regions of stability of paranematic (I), nematic 
(N) and log-rolling (L) states in the strain-rate-composition 
(a) and stress-composition (b) planes for L = 5.0. Note that 
the loop in (a) contains no stable states. Stability limits are 
calculated with respect to both order parameter and com- 
position fluctuations, for a given controlled stress. The thin 
loop in (b) encloses a region with no stable states, due to the 
instability of the composition equation. 



V. COMMON STRESS COEXISTENCE 



For common stress coexistence the interface lies in the 
velocity-vorticity plane, and inhomogeneities are in the 
y direction (see Fig. ||). The stress balance condition at 
the interface is <x • y uniform. a yy is taken care of by 
the pressure and a zy vanishes by symmetry (no flow in 
the z direction), leaving continuity of the shear stress 
a xy through the interface. The two coexisting phases / 
and // have strain rates and compositions partitioned 
according to 



4> = + (i - ()<t>u 
j = Cii + (i- Oin, 



(5.1) 
(5-2) 



where <j) and 7 are the mean composition and strain rate 
and £ is the fraction of material in phase /. 



A. Paranematic— flow aligning coexistence (l-N) 

Phase Diagram — Fig. || shows the tie lines computed 
on the {a xy — u) and (7— u) planes according to the pro- 
cedure outlined in Section III. Several features should 
be noted. Flow induces nematic behavior in what, in 
equilibrium, would be an isotropic phase. The tie lines 



are horizontal in the (a x 



plane, since phases coex- 



ist at a prescribed stress; and have a positive slope in 
the (7— it) plane because the more concentrated nematic 
phase flows faster. There is a critical point at sufficiently 
strong stress, whose existence is expected since the flow- 
aligning nematic and paranematic states have the same 
symmetry (Q is biaxial) and their major axes in the shear 
plane. 

More interesting is the changing slope of the tie lines. 
For weak stresses the equilibrium system is slightly per- 
turbed and the tie lines are almost horizontal. For high 
stresses the tie lines become more vertical and the com- 
position difference between the phases decreases. The 
slope of the tie lines determines the shape of the mean 
stress-strain-rate relation a xy (7) that would be measured 
in steady state experiments. 
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FIG. 9. Mean stress-strain-rate curves for coexistence at common stress, for L — 5.0 and A = 1.0. The solid lines denotes I 
and N branches; the dotted line in each figure denotes the stable N branch with which the I state coexists at the low strain rate 
boundary of the coexistence region, at a strain rate marked by an open circle o. The solid circles • and thick solid line denote 
the stress that would be measured in the banded regime. Phase coexistence occurs between phases of different compositions 
than the mean compositions (it = 2.555, 2.58, 2.685). The unstable portion of the homogeneous flow curve is shown in (a) and 
(b), but not (c). Note that the plateaus in the two- phase regions in (a) and (b) rather obviously do not satisfy an equal area 
construction with the underlying constitutive curve at the mean composition. 



Mean Constitutive Relations — Consider a composi- 
tion in the range where phase-separation occurs. For 
small applied stress 5 X y{i) varies smoothly until the two- 
phase region is reached. At this stress, a tiny band of 
high strain rate strongly-aligned nematic material ap- 
pears, with volume fraction determined by the lever rule, 
Eq. (5.1). The mea n strain rate 7 is determined by the 
lever rule, Eq. (5^), and the measured constitutive re- 
lation a X y(i) is non-analytic at this point (see Fig. |^). 
As the stress is increased further, the system traverses 
the two-phase region by jumping from tie line to tie line. 
Each successive tie line has a higher stress, a higher mean 
strain rate, and a steadily increasing volume fraction 
of nematic phase. The compositions of both coexisting 
phases change steadily through the two-phase region. 

The constitutive relation a xy (7) through the two phase 
region is determined by the spacing and splay of the tie 
lines. For mean compositions 4> close to the equilib- 
rium isotropic-nematic transition (Fig. ^c) the tie lines 
in the (7 — it) plane are fairly flat, so that the stress a xy 
changes significantly through the two-phase region; and 
the 'plateau' has definite curvature, reflecting the initial 
splay of the tie lines. For slightly lower mean composi- 
tions (Fig. ^|b) the 'plateau' is straighter and flatter, as 
can be seen in (Fig. [lO]), because the lines are more verti- 
cal in the (7— it) plane. Finally, for compositions near the 
critical point the plateau is flatter still but, more inter- 
estingly, phase coexistence occurs in a region where the 
stress-strain curve at the mean composition is no longer 
non-monotonic (Fig. [|a)! This is because stability in 
a two-phase system is also determined by the stability 
with respect to composition variations. In fact, the local 
chemical potential (J-(u) has negative slope and is unsta- 
ble on a segment of this curve. The tie line construction 
is a graphical expression of the explanation proposed by 
Schmitt et al. |2q], who attributed a sloped plateau to 



composition-dependence of the stress-strain constitutive 
relation. The general relation is given by 



da 



1 
VI 



1 



C 



JnVn JtVi 



(5.3) 



where m(a) is the slope of tie line with stress value a, the 
lines {a 1 (</)), on{4>), ji(<fi), jN{4>)} bound the phase coex- 
istence domains in the a—<fr and 7— <f> planes; rjk = dak/d'j 
is the local viscosity of the kth branch, and j' k = djk/dcj). 
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FIG. 10. a xy vs. 7 for common stress coexistence for 
L = 5.0 and A = 1.0. The solid lines connecting the high 
and low strain rate branches at each composition denote the 
composite flow behavior at coexistence. 

Measurements at controlled stress or controlled 
strain rate — Although these calculations are for phase 
separation at a common stress, one may perform experi- 
ments at either controlled stress or strain rate. All three 
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composite curves in Fig. ^ have similar shapes, so we ex- 
pect the same qualitative behavior for all compositions. 
Controlled strain rate experiments should follow the ho- 
mogenous flow curves, except for strain rates in the co- 
existence regime. Here we expect the steady state to 
eventually be the banded state. This should presumably 
occur by a nucleation event after some time, for start-up 
strain rates less than the I limit shown in Fig. ||a; and 
occur immediately for imposed strain rates beyond this 
stability limit. Conversely, upon decreasing the strain 
rate from the nematic phase we expect nucleated behav- 
ior for strain rates larger than the N limit, and insta- 
bility for smaller strain rates. In the metastable regime 
we expect the flow curve to follow the underlying ho- 
mogeneous constitutive curve for the given composition, 
until the nucleation event occurs. Interestingly, there is a 
small region (inside the loop in Fig. |^a) where the system 
is unstable when brought, at controlled strain rate, into 
this region from either the I or N states. This corresponds 
to constitutive curves with the multi-valued behavior of 
curve b in Fig. ||. 

Controlled stress experiments should exhibit similar 
behavior. Consider Fig. ^|b. For initial applied stresses 
larger than the minimum coexistence stress and less than 
the I limit of stability in Fig. |8|b, we expect the system 
to follow the homogenous flow curve until a nucleation 
event occurs. After nucleation the strain rate should in- 
crease, until either the proper plateau strain rate or the 
high strain rate nematic state is reached, depending on 
the magnitude of the stress. For stresses larger than the 
limit of stability we expect the system to become immedi- 
ately unstable to either a banded flow or a homogeneous 
nematic phase, depending on the magnitude of the stress. 
Metastability: Experiments — Experiments on worm- 
like micelles [p]Pj29| have found constitutive curves anal- 
ogous to those in, say, Fig. ||. In these experiments the 
plateau appears to be the stable states, while the por- 
tion of the constitutive curve (a 'spine') which extends 
to stresses above the onset of the stress plateau appears 
to be a metastable branch on which the system may re- 
main for a finite period of time under controlled stress 
or strain rate conditions. Refs. P, p9| , |5l| conducted con- 
trolled strain-rate experiments and found that the sys- 
tem follows the composite curves (without 'spines' that 
extend above the onset of the stress plateau) in Fig. ||, 
if care is taken to reach steady state. In these systems 
the plateaus were nearly flat, suggesting a very slight 
dependence of the flow behavior on composition. For 
controlled strain rate quenches into what corresponds to 
the two-phase region of Fig. ||a, the system took some 
time to develop shear bands and phase separate. This 
relaxation or 'nucleation' time decreased as the mean- 
strain rate was increased ^9|. It is not clear that they 
reached a limit of stability (which would be analogous to 
the I limit in Fig. ^ . The relaxation times were of order 
60 — 600 s, depending on temperature, mean composition, 
and mean strain rate. We emphasize that these experi- 
ments were on micellar solutions, which probably do not 



show an isotropic-nematic transition, but still display the 
same qualitative stress-strain-rate relationship as curve 
b in Fig. § 



Ref. [Q revealed different stress plateaus upon con- 
trolling either the strain rate or the shear stress (see 
Fig. 7 of Ref. [Q ) in cone-and-plate flow. In controlled 
stress experiments the stress plateau occurred at a stress 
of order 1.5 times the stress plateau observed under con- 
trolled strain rate conditions. Moreover, the flow curve 
under controlled stress conditions exhibited a stress max- 
imum and then a decrease in stress to an approximate flat 
plateau. One explanation for the high stress plateau un- 
der controlled stress conditions could be that the 'spine' 
never nucleated under controlled stress conditions, and 
the system smoothly transformed to the high strain rate 
phase. However, we do not have an explanation for the 
decrease and subsequent plateau in stress under applied 
strain rate conditions. 



In other experiments, controlled stress experiments 
revealed two kinds of metastable behavior [29). For 
a p < a < u jump, where a p is the minimum stress for the 
onset of banding in controlled strain rate experiments, 
the system maintained a strain rate on the 'metastable' 
branch for indefinite times (measured times were up to 
10 4 s). For a > a jump the system accelerated, after of 
order 10 3 s, and left the rheometer. For these systems 
it is not clear whether a stable high shear branch exists. 
An explanation for (Jj ump is lacking. Evidently the nu- 
cleation processes governing metastability at controlled 
stress and controlled strain rate are different. Clearly 
we need more experiments and theory about the nature 
of nucleation and metastability in controlled stress vs. 
controlled strain rate experiments. 



Polydispersity — Fig. |Tl| shows the effect of rod aspect 
ratio L on the phase diagram. A smaller rod aspect ra- 
tio couples more weakly to the flow, requiring a slightly 
larger strain rate to induce a transition to the nematic 
phase (Fig. [Tl]a) . The resulting stress is slightly smaller 
because, when the system enters the two phase region the 
stress is largely determined by that of the paranematic 
branch, which decreases with increasing L (Fig. [TT|b) . Al- 
though the equilibrium phase boundaries are close (see 
Eqs. ^9-3.10), the deviation is amplified considerably by 
applying flow. This suggests that flow enhances the nat- 
ural tendency of length polydispersity to widen biphasic 
regimes. 
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FIG. 11. Phase diagrams for L = 5.0 and L — 4.7 at 
common stress, for A = 1.0. 



B. Paranematic— Log rolling coexistence (l-L) 

Fig. [l2| shows the phase diagram calculated for coexis- 
tence between paranematic (I) and log rolling (L) states. 
As with l-N coexistence, the zero shear limit corresponds 
to the equilibrium biphasic region. However, for non- 
zero stress the biphasic region shifts in the direction of 
higher concentration. This is reasonable, since the sta- 
bility limit of the L phase shifts to higher concentrations 
with increasing stress (Fig. ^). Note also that, since the I 
and L phases have major axes of alignment in orthogonal 
directions, there is no critical point. Instead, the win- 
dow of phase coexistence ends when the I phase becomes 
unstable to the N phase. 

We have also computed phase coexistence between N 
and L phases. This occurs at much higher compositions 
(u > > 3.0) and has a narrow width in composition 
due to the very slight difference in viscosities of the two 
phases. Unfortunately, we cannot resolve this coexistence 
regime accurately within the numerical precision of our 
calculations and do not present these results here. 

The existence of two possible phase diagrams for com- 
mon stress phase separation raises an interesting ques- 
tion. Can one observe l-L coexistence? Notice that l-L 
coexistence can only occur for samples prepared at con- 
centrations at or above that necessary for equilibrium 
phase separation. One could prepare a phase separated 
isotropic-nematic mixture and, by wall preparation, field 



alignment, sedimentation, or other techniques, separate 
the phases into two macroscopic domains with the ne- 
matic phase in the log-rolling geometry. 
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FIG. 12. Phase diagram in the {o xy — u) and (7 — u) 
planes for paranematic-log-rolling coexistence, for L = 5.0 
and A = 1.0. The dotted lines are the limits of stability of the 
I and L phases (see Fig. |^) . 
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FIG. 13. Composite phase diagrams for l-L and l-N coex- 
istence at common stress for L = 5.0 and A = 1.0. We stress 
that this represents two overlayed phase diagrams, and not a 
single phase diagram. For example, there is no triple point 
implied by the intersection of the l-N and l-L phase diagrams. 

Upon applying shear, the system could then maintain 
coexistence and move through the l-L two-phase region. 
However, under controlled strain rate conditions, the I 
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material could decay into l-N coexistence (see Fig. 13). 
The resulting l-N coexistence occurs would quickly desta- 
bilize the entire l-L structure. Therefore the three-band 
structure N-l-L will not be present in this model, and it 
is probable that l-L coexistence could only exist under 
flow as a metastable state. Similar conclusions may be 
drawn by examining the phase diagrams in field-variable 
space, [i — a xy , as in Fig. |l8|a. In this case the chemical 
potential of the I phase, at l-L coexistence, is within the 
N region of the phase diagram for l-N coexistence, indi- 
cating a (possibly metastable) instability with respect to 
l-N phase separation. Moreover, the chemical potentials 
of the three phases are never the same, except at rest 
where the L and N states are identical apart from the 
rod orientations. 



VI. COMMON STRAIN RATE COEXISTENCE 



For coexistence at common strain rate the interface 
lies in the velocity-velocity-gradient plane, and inhomo- 
geneities are in the z direction (see Fig. 0). The stress 
balance condition at the interface is er • z uniform. As 
before, a zz is taken care of by the pressure while a yz and 
a xz are zero by symmetry (and because there are no sta- 
ble (fe components in the order parameter tensor) . With 
bands in the z direction, the strain rate in each band is 
set by the relative velocity of the two plates (or cylinders, 
in a Couette device), and the shear stresses differ. The 
mean applied stress a xy is the area average of the stress 
applied to each band. The coexisting phases have shear 
stresses and compositions partitioned according to 



4> = C^x + (1 - ()<t>n 

« y + (i - ck;, 



>xy 



(6.1) 
(6.2) 



where a xy is the mean shear stress. The interfacial equa- 
tions to solve are Eqs. (3. 27), (3. 30), and (3.31). 

Phase Diagram — Common strain rate l-N phase coex- 
istence is shown in Fig. |l4]. In this case the tie lines are 
horizontal in the (7— u) plane. They have a negative slope 
in the icr xy — u) plane because the paranematic I phase 
coexists with a denser and less viscous flow-aligning N 
phase. As with phase separation at common stress, there 
is a (very small) loop in the limits of stability in the con- 
trol variable plane (7— it) within which there are no stable 
homogeneous states. The careful reader will note that 
the limits of stability at a given stress (Fig. ||) are dif- 
ferent from the limits of stability at a given strain rate. 
This is phy sically correct, and will be discussed below in 

Sec. [vn q . 
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FIG. 14. Common strain rate phase diagram in the (7— u) 
(a) and (a xy — u) (b) planes, for L = 5.0 and A = 1.0. Also 
shown are the limits of stability of the I and N phases (calcu- 
lated for a given imposed strain rate, in contrast to Figures 

[| and |l3| in which the stability was calculated for an 
imposed stress.) 



There is an interesting crossover visible in the (o xy — u) 
plane. For higher mean compositions the fluid has a 
higher stress in its high strain rate one-phase region than 
in its low strain rate one-phase region; that is, respec- 
tively above and below the biphasic region in the Fig. 14 1. 
Conversely, for low enough compositions u < 2.67, the 
stress in the high strain rate region immediately outside 
the biphasic regime is actually less than the stress just 
before the system enters the biphasic region, as can be 
seen by the crossing of the solid and dashed phase bound- 
aries in Fig. |lj. 

This crossover is straightforward to understand. Since 
phase separation occurs at a given strain rate, and the 
stress of the N branch at a given composition and strain 
rate is always less than that of the corresponding I 
branch, we expect a decrease in the stress upon leav- 
ing the biphasic regime in cases where the coupling to 
composition is less important. We saw in the analysis at 
common stress that composition effects are less impor- 
tant (for l-N coexistence) at lower compositions and high 
strain rates, where the tie lines are more vertical. We ex- 
pect this near the critical point where the two phases be- 
come more and more similar. More generally, we expect 
this behavior in situations where phase separation occurs 
at a common strain rate into a shear-thinning state with 
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only slight changes in composition. In the more concen- 
trated regime, the coexistence plateau traverses a wider 
range of concentrations and strain rates, and emerges 
into the pure N phase with a higher stress (the width in 
strain rate of the phase coexistence regime is enough to 
overcome the shear thinning effect of the nematic phase). 
Mean Constitutive Relations — Figs. [l5|-^6| show the 
mean stress-strain-rate relations. As with common stress 
phase separation, the shape of the 'plateau' as the strain 
rate is swept through the two-phase region is not always 
fiat, and depends on the splay of the tie lines. At higher 
concentrations the plateau has a positive slope while, in 
accord with the crossover in the (a xy —u) phase diagram, 
for lower concentrations the plateau crosses over to neg- 
ative slope, which usually signifies a bulk instability. A 
simple argument, analogous to that for the stability of a 
bulk fluid, supports this. However, we note that a com- 



posite negative slope curve was accessed, and apparently 
found stable, by Hu et al. |(|| under controlled stress 
conditions. The negative slope in Fig. |l5| a is likely to be 
inaccessible under controlled stress conditions, and the 
instability argument may apply to controlled strain rate 
conditions. The general relation for the slope in the com- 
posite region is pq] 



da 



m C + vn (i - C) - m (i) 



vn(i - C) , Vi C 



'N 



(6.3) 



where miff) is the slope of the tie line with strain rate 7 
and a' k — dak /9<p. In the limit of no concentration dif- 

00), (7(7) is vertical through 



ference (S(j) = or 771(7) 
the two phase region. 
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FIG. 15. Mean stress-strain-rate curves for common strain rate coexistence for L = 5.0 and A = 1.0. The solid lines denote 
the stable I and N branches; the dotted line in each figure denotes the stable N branch with which the I state coexists at the 
low strain rate boundary of the coexistence region, at a strain rate marked by an open circle o. The solid circles • and thick 
solid line denote the stress that would be measured in the banded regime. The filled circles • and thick solid line denotes the 
stress measured under banded conditions. 



Measurements at controlled stress or controlled 
strain rate — For controlled strain rate measurements 
we expect behavior similar to that for phase separation 
at common stress. For start-up experiments with mean 
strain rates larger than the minimum strain rate for co- 
existence at a given composition, we expect the stress 
to follow the metastable branch until a nucleation event 
causes the stress to decrease to the plateau stress. The 
exception is a composition such as that in Fig. 15 1, for 
which the composite flow curve for l-N coexistence may 
be mechanically unstable. Similar results should apply 
upon decreasing the strain rate from the shear-induced N 
phase to below the N limit. As before, this expectation of 
a nucleation event is based on a possibly misguided anal- 
ogy with equilibrium systems which, nonetheless, is en- 
couraging given the experiments which see "nucleation" 
type behavior in micelles under flow 29 , 5l]] . 

For controlled stress, the situation is slightly different. 
For compositions with mean stress-strain-rate curves of 



the shape of Fig. |15|c, we expect similar behavior to that 
found for common stress phase separation. However, for 
compositions that yield curves such as Fig. |l5| a there is 
a window of stresses for which there are three possible 
states: homogeneous low strain rate and high strain rate 
branches, and a banded intermediate branch. We empha- 
size that we have not determined the absolute stability 
of any of these branches. A possibility is that the system 
has hysteretic behavior. For example, in start-up exper- 
iments the system would remain on the I branch until 
a certain stress, at which point it would nucleate after 
some time and transform to either the high strain rate 
N branch or coexistence. We cannot tell which state it 
might go to, from this analysis, but it seems likely that 
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it would jump straight to the N branch.''' The same be- 
havior (in reverse) would be expected upon reducing the 
stress from the high strain rate N phase. 
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FIG. 16. o xy vs. 7 for various compositions, for phase sep- 
aration at common strain rate and L = 5. 

Although there have been anecdotal reports of shear 
banding in the common strain rate geometry, there have 
been very few such results published. Bonn et al. ]3§| ] 
have recently reported results for sheared surfactant 
onion gels, along with visual confirmation of bands in the 
common strain rate geometry. In controlled strain rate 
experiments they found constitutive curves analogous to 
Fig. |l^a or |l5|b. In controlled stress experiments they 
found hysteretic behavior, with the system flipping be- 
tween high and low strain rate branches after some delay 
time, missing the coexistence 'plateau' regime. However, 
it is not clear that these were true steady state results. 

Stable 'negative-slope' behavior was seen in a shear- 
thickening systems which phase separates at common 
stress ]5^ , ]53t , under controlled stress conditions. In this 
case there was a single (mean) strain rate for a given ap- 
plied stress, and the measured constitutive relation had 
an S shape rather than the sideways S shape of Fig. FT|. 



VII. DISCUSSION 

A. Dependence on gradient terms 

Gradient terms appear in all equations of motion for 
K ^ and for any g, so to avoid unphysical equations 
without gradients (which cannot resolve interfaces) we 



must have A ~ g/K < 00. In the case of K = and fi- 
nite g the Q equation of motion has no explicit gradient 
terms and hence can, in principle, support discontinu- 
ous solutions. The <f> equation has gradients in this case, 
arisi ng f rom the term g (V</>) 2 in the free energy density, 
Eq. ( |3.2| ) , so the system will eventually reach a state with 
smooth solutions in both <f> and Q. Conversely, for g = 
there are gradient terms in both the Q and d> dynamics, 



with the latter arisi ng fr om the term 
energy density Eq. (3.2). 



1 (VQ) in the free 



Phase boundaries for a xy — 0.01,0.03 are shown in 
Fig. [l?]. For A £ (0.0 — 30.0) the phase boundaries are 
the same, within the precision of our numerical calcula- 
tions, while there is a distinct difference for A = 00. We 
have discretized the system on a mesh of 125 points, and 
the range of elastic constants is such that the width of 
the interface is at least 20 mesh points; large enough for 
smooth behavior and much smaller than the system size. 
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FIG. 17. I-N phase boundaries for common stress phase 
separation as a function of X/L, for L = 5.0. The diamonds 
♦ are for A = oo(K = 0,g = 1). 

We cannot rule out the possibility that changes in A 
shift the phase boundaries by small amounts below our 
accuracy, which is of order 0.1% in u, but the apparent 
independence of the phase boundaries on A is curious. 
One might be tempted to generalize and suggest that, for 
finite A, there exists a selection criterion which involves 
only the homogeneous equations of motion, rather than 
requiring the inhomogenous terms as in the interface con- 
struction. An interface construction may also be used to 
determine equilibrium phase boundaries, in which case a 
stationary interface is equivalent to minimizing a free en- 
ergy and the (relaxational) dynamical equations derived 



'If the system jumped from the I branch to the coextence 
branch, increasing the stress further would decrease the strain 
rate and return the system to the I branch. Since it originally 
nucleated from the I branch, it seems unlikely that the original 
jump could be to the coexisting plateau. 
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from a variational principle B . In the case of a van der 
Waals fluid this reproduces the Maxwell construction. 

A steady state equation for a single variable ip with 
homogeneous and inhomogencous terms of the form 



B. Which phase separation is preferred? 



Co = fhomW + finh(9ip/dy) 



(7.1) 



can be integrated to yield a solvability condition for o~q, 
which is equivalent to the stable interface method. In 
equilibrium /j n /, integrates exactly without an integrat- 
ing factor, since it typically arises from a variation of a 
free energy functional with respect to ip, and the result <7q 
(corresponding to the pressure in the van der Waals fluid) 
depends only on fhom- Out of equilibrium, integration is 
not so simple, and the solvability condition depends, gen- 
erally, on the form of the gradient terms |]60| , p5| . 

In the multivariable case considered here the steady 
state conditions for the order parameter and composi- 
tion are coupled differential equations which are not in- 
tegrable in shear flow. This is because of the terms 
K Q + Q k , t in Eq. ( fTl2| ) and (V 2 Q) Q - Q (V 2 Q) 
in Eq. ( p.30| ). In extensional flow k is symmetric, so that 
k -Q + Q ■ k t integrates to Tr(Q 2 /«)/2, while in shear 
flow this term can only be integrated by introducing an 
integral representation |7y] . Hence a first integral of the 
steady state equations cannot be found in shear flow, and 
it seems unlikely that a general condition involving only 
the homogeneous portion of the steady state equations 
can determine coexistence. While we appear to find, for 
this set of gradient terms, solvability conditions that are 
independent of A for A < oo, the relationship of this 
to a variational principle remains unknown. We have 
not exhausted the possible gradient terms. For exam- 
ple, higher order gradients in the free energy ((V 2 Q) 2 , 
etc.) would yield higher-order differential equations for 
the interfacial profile, and other square gradient terms 
such as Qq/jVq V/3<f> are possible [fnj. Hence, we believe 
that, for finite A, the apparent independence of our results 
on gradient terms only applies to the particular (simple) 
family of gradient terms we have chosen. The structure 
of the differential equations describing the steady states 
may change abruptly for A = oo, for which a term is lost 
in the differential equations, leading to a distinctly dif- 
ferent selection criteria and the shifted phase boundary 
in Fig. [It]. Unfortunately, this particular set of equations 
is too complex for this kind of analysis. For example, in 
a study of a simpler constitutive model, one can demon- 
strate that the selected stress depends on the detail form 
of the gradient terms |35| . 

Several workers have claimed to find an equal areas 
construction on the stress-strain-rate constitutive curve 
p0| . That is, the "plateau" as the system traverses the 
two-phase region is said to describe a path such that 
the areas above and below the plateau, enclosed by the 
plateau line and the underlying constitutive curve, are 
the same. This is not true here, as can be seen in Fig. UB. 



Having calculated both common strain rate and com- 
mon stress phase se par ation for the same system, and 
noticing from Figs. O and 14 that there are compo- 



sitions and shear conditions which lie inside the two- 
phase regions of all three calculated phase separations, 
we must address the question of which phase separation 
occurs. We have already argued that we expect l-L phase- 
separation at common stress to be metastable with re- 
spect to l-N phase separation at common stress. What 
about the relative stability of l-N phase separation at 
either common stress or common strain rate? 

With limited one-dimensonal calculations for systems 
of different symmetry (annular bands at common stress 
and stacked disk-like bands for common strain rate) it 
is impossible to calculate the stability of one interface 
profile with respect to another. Renardy calculated the 
stabilityof common stress coexistence to capillary fluctu- 
ations |p8| , which is a start; and such a stability analysis 
has been performed, in part, on the layer orientation of 
smectic systems in flow [ fl6| . However, some insight can 
be obtained by examining the "phase diagrams" in the 
chemical potential-field variable (either stress or strain 
rate) planes. The solid lines in Fig. [l8] are analogous to 
lines of phase coexistence in, for example, the pressure- 
temperature plane in a simple fluid. 



Consider Fig. [L8|a. Here, a xy and \x are the proper field 
variables for phase separation at a common stress, and 
the solid lines denote the l-L and l-N phase boundaries. 
The dashed line denotes the range of stresses at coexis- 
tence for common strain rate phase separation (Fig. 18 1), 
for which stress is a generalized density variable and 
strain rate the field variable. 



Fig. |l8| a indicates that, for a system undergoing l-N at 
a common strain rate, the chemical potential and stress 
for the I phase falls within the single phase N region of 
the common stress phase diagram. Hence, we expect this 
I phase to be unstable (or metastable) with respect to 
phase separation at common stress. Similarly, the con- 
trol parameters (chemical potential and stress) for the 
N phase coexisting at a common strain rate lie within 
the single phase I region for common stress phase separa- 
tion, which we also expect to be unstable (or metastable). 
Conversely, for a system coexisting at a common stress 
the I phase lies within the single phase I region of the com- 
mon strain rate phase diagram (Fig. p~8|j) , and similarly 
for the N phase. This suggests that phase separation at a 
common strain rate is unstable (or metastable) with re- 
spect to phase separation at common stress, while phase 
separation at a common stress is stable. 
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FIG. 18. Phase diagrams in the chemical potential /i vs. 
stress plane (a) and the /x-strain rate plane (b). The solid 
lines denote the phase boundaries for common stress phase 
coexistence in the /i— a plane (a) and for common strain rate 
coexistence in the fi — j plane (b). The dashed lines denote 
the coexisting stresses for common strain rate phase separa- 
tion within the common stress phase diagram (a); and vice 
versa in (b). 

Note that, ultimately, this selection of phase coex- 
istence geometries follows from the transition being a 
shear-thinning transition; for a shear thickening transi- 
tion the situation could be reversed. In this case phase 
coexistence at a common strain rate and a given fi would 
imply a shear-induced phase (analogous to the N phase) 
with a higher stress than the I phase. If the phase coex- 
istence line for common stress (strain rate) lay within 
a loop corresponding to the stresses (strain rates) for 
common strain rate (stress) coexistence, then common 
strain rate coexistence would be expected to be stable, by 
analogy with the isotropic-nematic shear thinning model. 
Obviously this argument is delicate. In a fluid where only 
one phase coexistence (either common stress or common 
strain rate) is supported by the dynamical equations this 
argument is moot. 
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FIG. 19. Phase diagrams in the (a) (i— 7 and (b) [i — a xy planes for I-N coexistence (the N state is stable for higher strain 
rate or stress, respectively). The thin vertical solid lines passing through B'-C and A-B-D denote phase coexistence at common 
strain rate and stress in (a) and (b), respectively. The broken lines marked 7 7 and iV 7 denote the coexisting states at common 
strain rate in the fi~a xy plane (b); while the broken lines I a and N a denote the coexisting states at common stress, in the 
y,— 7 plane, (c) is the mean stress vs. strain rate curve. Shown is a path A-B-C-D taken under the proposition that the system 
maintains a global minimum in chemical potential. Point A is at coexistence in the fJ.—<J xy plane (b), and hence corresponds to 
two points, on lines I a and N a , in the [i— 7 plane (a) for the two different strain rates of the coexisting phases. Similarly, point 
C corresponds to coexistence at common strain rate in (a), with the coexisting phases at different stresses lying on lines J 7 and 
N-y in (b), at the two points C. Points B and B' are coincident in (c), and correspond to a point switching from phase separation 
at common stress (B) to phase separation at common strain rate (B'). The path A-B-C in (c) may be traced in (b) by following 
the upper horizontal arrow until phase separation at common stress occurs at A, then along the segments A-B in (b) or A-B' 
in (a) until phase separation at common strain rate occurs at B ' . From this point until C the system phase separates along I-, 
and iV 7 , with a mean stress given by the thick diagonal solid arrow B-C in (b) and the thick segment B'-C in (a) . The system 
emerges from the two-phase region at C on N-y, and continues through D on the high strain rate branch. 



An alternative possibility is presented in Fig. |l9| if one 
argues that in steady state, among the possible phases 
which are compatible with the interface solvability con- 
dition, the chemical potential reaches its absolute min- 
imum. Consider increasing the strain rate for a given 
mean concentration. The thick horizontal arrows in 
Fig. [l9|a-c denote the [i(a X y) and /i(7) paths for the ho- 
mogeneous high and low shear rate states, in the two 
phase diagrams. In Fig. |l9| a the path is A-B'-C-D, in 
Fig. jlf^ b the path is A-B-C-D, and in Fig. [l9| c the path is 
A-B-C-D. 

This path ensures that the system maintains the mini- 
mum chemical potential for an imposed strain rate. Upon 
increasing the strain rate from zero the system remains 
in the one phase region until A is reached, at which point 
phase separation at common stress occurs. Note that 
A spans two points of coexistence in the \x — 7 plane 
(Fig. |l9|a) on the lines I a and N a . Upon further increas- 
ing the strain rate the system continues to phase separate 
at common stress, following the segment A-B in Fig. p~9| b 
and the two (coexisting) segments A-B in Fig. [l9|a. The 
mean chemical potential and strain rate follow the di- 
agonal segment A-B' in Fig. |l9|a. Upon increasing the 
strain rate above B, the system can continue to maintain 
its lowest chemical potential by phase separating with a 
common strain rate in the two phases, i.e. with shear 
bands in the vorticity direction. Hence, the system next 
follows the path B ' -C in Fig. |l9a f /x-^y plane) and the two 
coexisting paths B'-C in Fig. |1 9|z> (fi — cr X y plane), with 
the mean chemical potential and stress following the di- 
agonal segment B-C in Fig. [l9|b. Finally, upon increasing 
the strain rate above C the system continues along the 
high strain rate branch. The thin diagonal lines with ar- 
rows, B ' -D in Fig. |l9| a and B-D in Fig. |l9| b show the path 
that would be taken if the system passed through the 
two-phase region entirely with a common stress in the 
two phases. This scenario follows from minimizing the 
chemical potential, and its correctness, of course, should 
be further examined by the full time evolution of the 
original dynamic equations. 

It is probable that boundary conditions also play a role. 
Consider a Couette device. Typically the walls provide 
uniform boundary conditions in the azimuthal direction, 



while the slight inhomogeneity of Couette flow induces an 
asymmetry between the inner and outer cylinders. The 
slightly higher stress near the inner wall provides a pref- 
erence for the high strain rate nematic phase, and hence 
might enhance the stability of common stress phase sep- 
aration. Similarly, the intrinsic inhomogeneity (although 
weaker) in cone-and-plate rheometry induces a preference 
for the common stress interfacial configuration ]46| ] . 

We are also unable to say anything about the num- 
ber or spacing of bands. Analogies with equilibrium sys- 
tems suggest that phase separation would coarsen until 
the system formed two bands at different strain rates 
(for phase separation at a common stress). This is the 
behavior seen in visualizations of flow in Couette, cone- 
and-plate, and pipe geometries, where the intrinsic in- 
homogeneity provides a "seed" for macroscopic phase 
separation @,E3 [|^^2| • Recent visualization of banding 
in lamellar surfactant systems p8[ indicate that phase 
separation at a common strain rate can exhibit bands 
(disc-like bands in Couette flow) whose initial spacing 
depends on the applied strain rate and coarsen in time. 
Unlike common stress bands, which are expected to (and 
do) form macroscopic bands in Couette flow, there is no 
boundary effect in Couette (aside from perhaps sedimen- 
tation) which would encourage common strain rate bands 
to coalesce readily. Normal stresses may play a role in 
this process. 



C. Stability at prescribed stress or prescribed strain 

rate 



In calculating the phase diagrams we have calculated 
the stability of the fluid under conditions of either fixed 
strain rate or fixed stress. These limits of stability, anal- 
ogous to spinodals in equilibrium systems, are displayed 
in Figures || and [l4|. Note that the stability limits and 
critical points differ, depending on the control variable 
(stress a or strain rate 7). To see why this is, note that 
schematically the dynamical equations of motion have 
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the form 



djx = f(x, 7) 



(7.2) 
(7.3) 



where x comprises the dynamical variables (order param- 
eter Q and composition <f)). The second equation relates 
the stress to strain rate and dynamical variables at steady 
state (or in the zero Reynolds number limit), which im- 
plies that the strain rate 7, for a given stress, is a function 
7 = 7(0-, x). Consider fluctuations about a steady state 
xo: x = <5x + xq. The dynamics for the fluctuation obeys 

= {M 7 + 6M a } ■ 5x (7.5) 

The limit of stability for common strain rate is calcu- 
lated using the fluctuation matrix M 7 , while the limit of 
stability for common stress was calculated using M 7 + 
<5M CT . These correspond to different stability criteria. 

The question of which spinodal could be observed in 
an experiment relies on the accuracy of prescribed stress 
and prescribed strain rate rhcomcters. For a rheometer 
operating at a prescribed s train rate, then if (foe goes un- 
stable t hrou gh M 7 (in Eq. 7.2), the stress increases due 
to Eq. (7.3) and no attempt is made to control it, lead- 
ing to instability. However, consider a rheometer which 
mainta ins a prescribed stress. If the system goes uns tabl e 
in Eq. (7.2), the bulk stress will change due to Eq. (7.3). 
A sensitive and fast enough rheometer will respond by 
adjusting the strain rate accordingly, to maintain the im- 
posed stress. Hence, instability would be determined by 



the sum M, 



6M r 



Similarly, in equilibrium systems a locus of stability 
may be defined by, for example, the diverging of isother- 
mal or adiabatic, or isobaric or isochoric, response func- 
tions (or the vanishing of the appropriate modulus). For 
example, the isothermal and adiabatic compressibilities 
Kt and K$ differ by a term proportional to the quotient 
of the square of the thermal expansion coefficient a p and 
the isobaric heat capacity c p : 



Kn 



K S = 



vTal 



(7.6) 



where v is the specific volume. K T X vanishes along the 
spinodal line v s (T), while it is evident that Kg 1 (propor- 
tional to the sound speed) does not. However, in equilib- 
rium, the critical point is uniquely defined in phase space, 
which is related to the fact that, for example, pressure 
is a unique function of the volume, and is in fact a state 
variable. Conversely, we can see from the shape of the 
stress-strain rate curves for the Doi model (e.g. Fig. ||c- 
e) , that the stress can be a multivalued function of strain 
rate; i.e. it is not a state function. Hence there is no 
compelling reason to expect critical points at imposed 
strain rate to be the same as critical points at imposed 



stress. Similarly, the true spinodal, or locus of stabil- 
ity, is uniquely defined in an equilibrium system because 
of the convexity requirement on the entropy |72] ], and 
there is no such universal convexity requirement (bar- 
ring entropy production, which is minimized only under 
restricted conditions, and only locally rather than glob- 
ally) for non-equilibrium systems. 



D. An analogy with equilibrium systems? 

The liquid crystalline suspension under flow, indeed 
any system which undergoes a macroscopic bulk flow- 
induced phase transition, is analogous to an equilibrium 
ternary system comprising species A, B, and solvent. In 
our case, the roles of A and B are played by the rigid 
rod composition (f> and either the stress a or strain rate 
7, depending on the nature of the phase separation. For 
phase separation at common stress, the phase diagram 
in the stress-composition plane a — (f> is analogous to the 
fJ-A — <Pb plane for the equilibrium system, while the 7— </> 
plane is analogous to the </>a — 4>b plane. In either case, 
the density variables, {7, <fi} in flow and {4>a, 4>b} in the 
analogous equilibrium system, are different in the two co- 
existing phases. The slope of the "plateau" in the a — 7 
plane, as the system traverses the two-phase region of 
the phase diagram, is analogous to a slope in the fiA—<t>B 
plane, the latter indicating that the chemical potential 
(or osmotic pressure) of the two phases varies across the 
coexistence region. 

Can this analogy be extended to the possibility of 
phase separation at common stress or common strain 
rate? Certainly, one can consider a ternary system under 
conditions of either imposed <Pa or imposed \xa , for which 
one generally expects difference spinodal lines. That is, 
the spinodal is determined by the instability of a matrix 
in the two-dimensional space spanned by 4>a and </>b, 
and fixing 4>a or [ia projects this instability onto differ- 
ent subspaces. Experimental conditions may dictate that 
the spinodal line under fixed <Pa is more likely to be seen 
by, since 4>a is conserved and cannot equilibrate quickly 
to satisfy an imposed \x a- However, we are not aware of 
any ternary equilibrium system for which the equilibrium 
coexistence conditions can differ; that is, equilibrium is 
always specified by equality of fiA and /j^ in the two 
phases, and never by equality of 4>a- 



VIII. SUMMARY 

In this work we have proposed a straightforward phe- 
nomenological extension to the Doi model for a solution 
of rigid rod particles. We have added entropic terms, and 
included inhomogeneous terms in order to calculate, for 
the first time, phase separation in shear flow. The main 
results of this study are as follows: 
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1. Phase separation may occur under conditions of 
common stress or common strain rate, with differ- 
ent interface orientations with respect to flow ge- 
ometry for the two cases. 

2. Although both phase separations are possible, the 
phase diagrams in the /i — a xy and fi — 7 planes 
(Fig. |l8|) suggest that phase separation at a com- 
mon strain rate is metastable. This can be traced, 
for this model, to the shear-thinning character of 
the transition. For a shear thickening transition 
an equivalent argument suggests that (if both are 
kinematic possibilities) common stress phase sep- 
aration is metastable with respect to strain rate 
phase separation. 

3. The limits of stability ( "spinodals" ) and critical 
points for systems at prescribed stress and pre- 
scribed strain rate differ; the difference of spinodals 
is similar to equilibrium behavior, while the differ- 
ence of critical points is related to the fact that 
neither stress or strain rate are always unique state 
functions. 

4. An argument based on minimizing the chemical po- 
tential predicts a complex crossover from common 
stress to common strain rate phase separation for 
controlled strain rate experiments. The veracity of 
this assumption is unknown. 

5. We have calculated phase coexistence among three 
phases (paranematic I , flow-aligning nematic N , and 
log-rolling nematic L), where only two phases ex- 
isted in equilibrium. We expect l-N phase coex- 
istence to be the stable configuration (Fig. 13), 
although l-L phase coexistence could exist as a 
metastable state with approprate preparation. We 
do not expect three phase coexistence for this 
model. 

6. We have demonstrated how to calculate the mean 
stress-strain-rate relationship a(j) in the coexis- 
tence region. The shape of 17(7) is determined by 
the composition and strain rates of the coexisting 
phases p3]. 

7. A phase-separated system can exhibit an appar- 
ently unstable constitutive relation, with negative 
slope daxy/dj. Experiments have accessed such 
negative slope composite curves under controlled 
stress (rather than controlled strain rate) condi- 
tions @. 

8. Our method of solution is general and relies on the 
existence of a set of dynamical equations of motion 
for the structural order parameter of the particu- 
lar transition, including the dynamic response to 
inhomogeneities . 

9. For A = g/K finite the phase boundaries we have 
found are, within our accuracy, independent of the 



relative magnitude of the gradient terms in our 
free energy. Although this suggests that, for the 
restricted set of inhomogeneities we have incorpo- 
rated, a selection criterion exists involving only the 
homogeneous equations of motion, this is not true 
in general for complex fluids in flow | j3o| . For A = 00 
the phase boundaries are slightly shifted. 

10. Studies at different aspect ratios suggest that shear 
flow enhances polydispersity effects relative to their 
effect on equilibrium phase boundaries. 

We close by enumerating several open questions. First, 
systems such as wormlike micelles probably possess some 
combination of a perturbed isotropic-nematic transition 
and a dynamic instability of the molecular constitu- 
tive relation. It is conceivable that suitable composi- 
tions of these systems could yield a stress-strain-rate- 
composition surface (Fig. ||) with multiple folds. Second, 
it would be desirable to have a model shear-thickening 
system in which to calculate properties of banded flows, 
to compare and contrast with the shear-thinning system 
studied here and to understand experiments on a wide 
range of systems, including clays and surfactant systems. 
Third, we have not addressed the number and possible 
coarsening of bands and band configurations; and the 
kinetics of phase separation has hardly been treated the- 
oretically with experiments also at an early stage 
p|, p9|j5l| . Finally, we do not yet know the conditions 
which may, if at all, distinguish between common stress 
or common strain rate phase coexistence. 
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